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We examine indirect signals produced by neutralino self-annihilations, in the galactic halo or 
inside celestial bodies, in the frame of an effective MSSM model without gaugino-mass unification 
at a grand unification scale. We compare our theoretical predictions with current experimental 
data of gamma-rays and antiprotons in space and of upgoing muons at neutrino telescopes. Results 
are presented for a wide range of the neutralino mass, though our discussions are focused on light 
neutralinos. We find that only the antiproton signal is potentially able to set constraints on very 
low-mass neutralinos, below 20 GeV. The gamma-ray signal, both from the galactic center and 
from high galactic latitudes, requires significantly steep profiles or substantial dumpiness in order 
to reach detectable levels. The up-going muon signal is largely below experimental sensitivities for 
the neutrino fiux coming from the Sun; for the fiux from the Earth an improvement of about one 
order of magnitude in experimental sensitivities (with a low energy threshold) can make accessible 
neutralino masses close to O, Si and A4g nuclei masses, for which resonant capture is operative. 

PACS numbers: 95.35.+d,98.35.Gi,11.30.Pb,12.60. Jv,95.30.Cq 



I. INTRODUCTION 

In Supersymmetric models without gaugino-mass uni- 
fication at the grand unification scale, neutralinos can be 
lighter than the current lower bound of 50 GeV, which 
instead occurs in the case of gaugino-universal models. 
In Refs. 0, we discussed the properties of these light 
neutralinos as relic particles (i?-parity conservation is 
assumed) and showed that an absolute lower limit of 6 
GeV on the neutralino mass to^ can be placed by apply- 
ing the most recent determinations of the upper bound on 
the Cold Dark Matter (CDM) content in the Universe, in 
combination with constraints imposed on the Higgs and 
supersymmetric parameters by measurements at colliders 
and other precision experiments, like the muon anoma- 
lous magnetic moment and the rare decay 6 — > s -I- 7. In 
Rcfs. Q 01 we also showed that direct detection rates 
for light relic neutralinos make these particles detectable 
with WIMP direct search experiments with current tech- 
nologies. A comparison of our predictions with interven- 
ing experimental results was presented in Ref. Q . 

In the present paper we examine light neutralinos in 
connection with the indirect signals which can be pro- 
duced by neutralino self-annihilations in the galactic halo 
or inside celestial bodies. We compare our theoretical 
predictions with current experimental data on measure- 



ments of gamma-rays and antiprotons in space and of 
up-going muons at neutrino telescopes. Results are pre- 
sented for a wide range of the neutralino mass, from the 
established lower bound of 6 GeV up to 500 GeV. How- 
ever, our discussions are focused on light neutralinos (i.e. 
neutralinos with < 50 GeV), since these are not usu- 
ally considered in current literature 0. 

The structure of the paper is the following. In Sec. ^ 
we briefly summarize the gaugino non-universal super- 
symmetric model and the properties of light neutralinos 
which arise in this framework. In Sec. IIIII we discuss 
the dark matter density distribution in the galactic halo, 
which is relevant to indirect detection signals, especially 
to the gamma-ray flux. In Sec. II VI we present the calcu- 
lation and comparison with data of the gamma-ray flux: 
we consider both the signal coming from the galactic cen- 
ter and from high galactic latitudes. In Sec. I^we discuss 
the antiproton signal, whereas in Sec. I VII we show our re- 
sults for the indirect signals at neutrino telescopes. Our 
conclusions are drawn in Sec. IVIII 



II. SUPERSYMMETRIC MODELS WITHOUT 
GAUGINO-MASS UNIFICATION 



A typical assumption of supersymmetric models is the 
unification condition for the three gaugino masses Mi^2,3 
at the GUT scale: Mi = M2 = M3. This hypothe- 
sis implies that at the electroweak scale Mi ~ 0.5 M2. 

♦Elec tronic address: [ bottinoii^to.infn.it, donatol^to.inl'n.it, fornengoti toSmW, (iiiAifclwiS'tifti.W fl Condition the bound on the ncu- 
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is derived from the lower bound on the chargino mass 
(which depends on M2 but not on Mi) determined at 
LEP2: m^± > 100 GeV. By allowing a deviation from 
gaugino-universality, the neutralino can be lighter than 
in the gaugino-universal models when Mi = R M2, with 
R < 0.5. In this case current data from accelerators do 
not set an absolute lower bound on m^. 

We consider here an extension of the MSSM which 
allows for a deviation from gaugino-mass universality by 
the introduction of the parameter R, varied here in the 
interval: (0.01 0.5). This range for R implies that the 
accelerator lower bound on the neutralino mass can be 
moved down to few GeV for R ~ 0.01. The ensuing light 
neutralinos have a dominant bino component; a deviation 
from a pure bino composition is mainly due to a mixture 
of B with Hi 0,0,1^. Notice that our range of R includes 
also the usual model with gaugino-mass universality. 

We therefore employ an effective MSSM scheme at the 
electroweak scale, defined in terms of a minimal num- 
ber of parameters, only those necessary to shape the 
essentials of the theoretical structure of MSSM and of 
its particle content, supplemented by the gaugino non- 
universality parameter R. The assumptions that we im- 
pose at the electroweak scale are: a) all squark soft- 
mass parameters are degenerate: rriq. = niq] b) all slep- 
ton soft-mass parameters are degenerate: mj = mf 
c) all trilinear parameters vanish except those of the 
third family, which are defined in terms of a common 
dimensionless parameter A: Aij^ = = Arriq and 
Af = Amj. As a consequence, the supersymmetric pa- 
rameter space consists of the following independent pa- 
rameters: M2, /i, tan (5, rriA, nig, nij, A and R. In the pre- 
vious list of parameters we have denoted by fi the Higgs 
mixing mass parameter, by tan /3 the ratio of the two 
Higgs v.e.v.'s and by niA the mass of the CP-odd neutral 
Higgs boson. 

In the numerical random scanning of the supersymmet- 
ric parameter space we have used the following ranges: 
1 < tan/3 < 50, 100 GeV < M2 < 1000 GeV, 
100 GeV < mq^mj < 1000 GeV, sign(/z) = -1,1, 
90 GeV <mA< 1000 GeV, -3 < A < 3, in addition to 
the above mentioned range 0.01 < R < 0.5. We impose 
the experimental constraints: accelerators data on super- 
symmetric and Higgs boson searches and on the invisible 
width of the Z boson, measurements of the branching ra- 
tio of the 6 s + 7 decay and of the upper bound on the 
branching ratio of Bg — > /i+ + fi^ , and measurements of 
the muon anomalous magnetic moment = [g^^ — 2)/2. 
The range used here for the & ^ s + 7 branching ra- 
tio is 2.18 • lO"'' < BR{b s + 7) < 4.28 • 10"^ 
0. For the branching ratio of i?s ^ ^+ + jjr we em- 
ploy the upper limit BR{Bs ^ + n^) < 7.5 • 10""^ 
(95% C.L.) 0; for the theoretical evaluation we have 



used the results of |3 with inclusion of the QCD radia- 
tive corrections to the b Yukawa coupling jg]. For the 
deviation of the current experimental world average of 
flp from the theoretical evaluation within the Standard 
Model we use the 2a range: -142 < Aa^ • 10" < 474; 
this interval takes into account the recent evaluations of 
Refs. 1^ We notice that gluinos do not enter di- 
rectly into our loop contributions to BR{b — > s -|- 7) and 
BR{Bs — > /i^ -l- A*")! since we assume flavor-diagonal 
sfermion mass matrices. Thus, gluinos appear only in 
the QCD radiative corrections to the b Yukawa coupling; 
in the evaluation of these effects the value of the rele- 
vant mass parameter M3 is taken at the standard unifi- 
cation value M3 = M2 az{Mz)/a2{Mz), where a^lMz) 
and a2{Mz) are the SU(3) and SU(2) coupling constants 
evaluated at the scale Mz- 

The new data on the cosmic microwave background 
from WMAp[ll|. used in combination with other cosmo- 
logical observations, mainly galaxy surveys and Lyman-a 
forest data, are sharpening our knowledge of the cosmo- 
logical parameters, and in particular of the amount of 
dark matter in the Universe. From the analysis of Rcf. 
\n\ , we obtain a restricted range for the relic density of a 
cold species like the neutralinos. The density parameter 
of cold dark matter is bounded at 2a level by the values: 
{^CDMh'%,in = 0.095 and {flcDMh^Ue.. = 0.131. This 
is the range for CDM that we consider in the present 
paper. An independent determination for the content 
of cold dark matter in the Universe is provided by the 
Sloan Digital Sky Survey Collaboration this new 

data agrees with the results of Ref . • 

We recall that the relic abundance Jl^^/i^ is essentially 
given by fl^h"^ oc< a^^nnV >^Ji, where < a^^nnV >int is 
the thermally-averaged neutralino annihilation cross sec- 
tion times average velocity, integrated from the freeze- 
out temperature Tf to the present one Tq. The quantity 
CTann cnters also in the calculation of the indirect signals 
that will be discussed in the following sections. In the 
evaluation of Cann we have considered the full set of avail- 
able final states: fermion-antifermion pairs, giuon pairs, 
pairs of charged Higgs bosons, one Hig gs boson and one 
gauge boson, pairs of gauge bosons [l3j. We have not in- 
cluded coannihilation in our evaluation of the neutralino 
relic abundance, since, at variance from a constrained 
SUGRA scheme, in our effective supersymmetric model 
a matching of the neutralino mass with other masses is 
accidental, i.e. not induced by some intrinsic relationship 
among different parameters. Introducing coannihilation 
would only produce an insignificant reshuffle in the rep- 
resentative points of the scatter plots displayed in the 
present paper, without a modification of their borders, 
which are the only feature of physical significance. 

The relic abundance fl-^h'^ of neutralinos lighter than 
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50 GeV which arise in our class of gaugino non-universal 
models has a relatively simple structure in terms of dom- 
inant diagrams in the annihilation cross section P, 
Here we just remind that combining our calculation of 
the relic abundance of light ncutralinos with the value 
of (i^cDA/^^)maxi absolutc lower bound on the ncu- 
tralino mass of 6 GeV can be set Q, Q- We note 
that within our present scanning of the supersymmet- 
ric parameter space, the lower limit on the ncutralino 
mass shifts to about 7 GeV, when the upper bound on 
BR{Bs < 7.5-10"^ is implemented (this con- 

straint was not included in It is remarkable that 

a lower limit on is set not by searches at accelerators, 
but instead by cosmological arguments. 



III. DARK MATTER IN THE GALAXY 

Signals due to ncutralino self-annihilation in the halo 
depend quadratically on the dark matter density distri- 
bution p(r), and are therefore very sensitive to the fea- 
tures of this physical quantity. Two properties are of 
special relevance: 1) the behavior of the density distri- 
bution in the galactic center (GC), 2) the extent of the 
density contrast (dumpiness), which represents the de- 
viation of the actual density distribution from a smooth 
distribution. 

The most commonly used density distributions can be 
parametrized by the following spherically-averaged den- 
sity profile 01 : 



p(r) 



1 + [r^/ar 



(/3-7)/q 



(r/a)" J ' 

where r ~ \ f\, tq = 8 kpc 0| is the distance of the Sun 
from the galactic center, a is a scale length and pi is the 
total local (solar neighborhood) dark matter density. In 
particular, the isothermal density profile corresponds to 
(a, (3, 7) = (2,2,0), the Navarro, Frenk and White (NFW) 
profile 0] corresponds to (a,/?, 7) = (1,3,1) and the 
Moore et al. profile [l3 to (a,/3,7) = (1.5,3,1.5). The 
two latter profiles, both derived from numerical simula- 
tions of structure formation, differ noticeably in their be- 



havior at small distances from the GC: 



for the NFW 



and r~^'^ for the Moore et al. profile, with ensuing large 
differences in the size of the expected signals for WIMP 
annihilation from the central region of the Galaxy. 

Recent results of extensive numerical simulations, 
aimed at an analysis of the inner structure of KCDM 
halos. strongly disfavor a behavior as singular as r~^'^, 
but also indicate that a NFW profile is likely not to be 
adequate at small distances from the GC hj. It turns 
out that the density profile is not described by a singular 
power-law at small distances; rather, in this asymptotic 



regime, the numerical results are well fitted by a profile 
whose logarithmic slope 5{r) = — (i(ln p(r))/d(ln r) is 
given by 5{r) ^ r^" with a ~ 0.17. This leads to a non- 
singular dark matter density distribution function of the 
form [13: 



p{r) 



p-2 CXp 



r 
r-2 



1 



(2) 



where r_2 is the radius where the logarithmic slope is 
6 ~ —2. and p_2 = p(?'-2)- These various distributions 
mainly differ in their behavior at small r. We wish to 
stress here that anyway current cosmological simulations 
are not reliable for radii smaller than an r^i^ ~ 0.1 - 
1 kpc. We also notice that singular profiles are subject 
of debate in current literature, with analyses pointing 
to inconsistencies with observational data on rotational 
curves [lij . 

Furthermore we recall that other density profiles are 
able to describe the dark matter halo, including for in- 
stance different classes of logarithmic and power-law po- 
tentials, axisymmetric distributions or even triaxial dis- 
tributions j20j. In the following we will concentrate on 
the standard isotropic density profiles of Eq. (Q) and on 
the new profile of Eq. (|2Jl deduced from numerical sim- 
ulations. For definiteness, we use as a reference model 
the NFW density distribution, and discuss the deviation 
from this reference case when the other density profiles 
are considered. Most signals do not depend or are only 
mildly dependent on the critical behavior of the dark 
matter density around the GC: this occurs for the neu- 
trino fluxes from the Earth and the Sun, for antiprotons 
and for gamma-rays coming from large galactic latitudes. 
On the contrary, gamma-rays coming from the GC are 
very sensitive to the inner parts of the Galaxy and the 
differences between the different halo profiles will be ex- 
plicitly discussed. 

A key parameter for all the density distributions is the 
value for the total local dark matter density pi. This 
parameter can be determined for each density profile 
assuming compatibility with the measurements of rota- 
tional curves and the total mass of the Galaxy l^^. For 
instance, a simple modelling of the visible and dark com- 
ponents of the Galaxy showed that pi can range from 0.18 
GeV cm"'^ to 0.71 GeV cm^'^ for an isothermal sphere 
profile, from 0.20 GeV cm~'^ to 1.11 GeV cm~'^ for a 
NFW distribution and 0.22 GeV cm'^ to 0.98 GeV cm-^ 
for a Moore et al. shape [i^. For definiteness, our re- 
sults will be presented for pi = 0.3 GeV cm"'^ for all the 
density profiles employed in the present analysis. The 
parameter pi enters as a mere scaling factor in the signal 
fluxes: the effect of varying pi is therefore easily taken 
into account. 

Once the density profile which describes the total dark 
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matter density in the galactic halo is chosen, the actual 
neutralino density distribution is taken to be: 



(3) 



where ^ accounts for the fact that neutralinos could be 
only a fraction of the total cold dark matter < 1). This 
characteristic is linked to the actual relic abundance of 
neutralinos and is accounted for by using the standard 
rescaling prescription: ^ = miii[l^ flv/i^/(f2cr>M^^)min]- 

As was noticed in Refs. [2l[|22,|23 an effect of density 
contrast in the dark matter distribution could produce a 
strong enhancement effect in signals due to x~X annihi- 
lations in the halo. This property was subsequently con- 
sidered in connection with various signals (gamma rays, 
positrons, antiprotons) sometimes under the as- 

sumption of a strong dumpiness effect, at the level of a 
few orders of magnitude. However, according to a recent 
analytical investigation on the production of small-scale 
dark matter clumps , the dumpiness effect would not 
be large, with the result that the ensuing enhancement 
on the annihilation signals is limited to a factor of a few. 
Similar conclusions are also reached in Ref. by using 
results of high-resolution numerical simulations. 



IV. GAMMA RAYS 

The flux of gamma rays originated from 

neutralino pair annihilation in the galactic halo j2ll . E3 . 
and coming from the angular direction ip is given by: 



1 (cTannW) dN^ 1 



dE. 



(4) 



where ((Tannic) is the annihilaton cross section times the 
relative velocity mediated over the galactic velocity dis- 
tribution function and dN^/dE^ is the energy spectrum 
of 7-rays originated from a single neutralino pair annihi- 
lation. The quantity /(V') is the integral of the squared 
dark matter density distribution performed along the line 
of sight: 



(5) 



1.0. s 



and -0 is the angle between the line of sight and the 
line pointing toward the GC. The angle ip is related to 
the galactic longitude I and latitude h by the expres- 
sion cos ip = cos I cos b. A point at a distance A from 
us and observed under an angle %p is therefore located at 
the galactoccntric distance r = y^A^ + Tq — 2Aro cos -0. 
The factor of 1/2 in Eq.|@} is due to the fact that the 
gamma-ray flux depends on the number of neutralino 



pairs present in the galactic halo, as pointed out in Ref. 
|29j. This factor of 1/2 applies as well to any other indi- 
rect detection signal which depends on the annihilation 
of a pair of Majorana fermions in the galactic halo, like 
positrons, antiprotons, antideuterium. In the case of dark 
matter composed of Dirac fermions, the statistical factor 
would instead be 1/4, if p{r) describes the total dark 
matter distribution ascribed to the given Dirac species. 



A. The source spectrum 

As far as the annihilation of light neutralinos is con- 
cerned (namely for neutralino masses below the thresh- 
olds for gauge-bosons, higgs-bosons and t quark produc- 
tion), the production of 7-rays in the continuum takes 
contribution mainly from the hadronization of quarks 
and gluon pairs produced in the neutralino annihila- 
tion process. The subsequent ir^ production and decay 

TT*^ > 27 give usually the dominant contribution. In 

this case the 7-ray energy spectrum is given by: 



dN^ 
dK 



PiE^,E^)^dE^, (6) 



where P{E^,E^) = 2{El - to^)-^/^ is the probability 
per unit energy to produce a 7-ray with energy E^ out 
of a pion with energy E"^, while dN^^/dE^^ is the pion 
yield per annihilation event. 

We have evaluated the quantity dN^^/dE.,,- by means 
of a Monte Carlo simulation with the PYTHIA pack- 
age [s^- The Monte Carlo code has been run by in- 
jecting qq and gluon pairs back-to-back at fixed ccnter- 
of-mass energy Ecm = 2 . Since quarks and gluons 
are confined, they contribute to a complex final-state 
pattern of out-going hadronic strings decaying to phys- 
ical hadrons through fragmentation. In the Lund string 
scheme, fragmentation is an intrinsically scale invariant 
process. This implies that, in the rest frame of the de- 
caying string, the final-state spectrum is invariant in the 
variable x = Ep /mstring, where Ep is the energy of the 
given final state and mgtring is the total string mass. Were 
it not for showered gluons, a qq pair from neutralino an- 
nihilation would produce, in the reference frame of the 
two annihilating neutralinos, a single hadronic string at 
rest with rngtring = 2to;^, subsequently fragmenting to 
produce (among other particles) a pion spectrum which 
would be scale-invariant in the variable y = E-^/m^. 
However, due to showering, this scale invariance is signif- 
icantly broken, since the pion energy-spectrum is given 
by the superposition of different decaying strings boosted 
at different energies. Therefore the pion spectrum at a 
given neutralino mass cannot be obtained from that cal- 
culated at a different rriy. 
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We have therefore evaluated the pion yield per anni- 
hilation event, dN^ {ni-^, E^^) / (IEt^ , for each final state 
F = ff,gg, at different neutralino masses: = 6, 10, 
50, 100, 500 and 1000 GeV and for pion energies ranging 
from ttIt^o to in 100 equal bins in logarithmic scale. 
In order to optimize our numerical calculations, we then 
obtain the pion spectrum at neutralino masses and pion 
energies different from the ones sampled through a two 
dimensional numerical interpolation. Wc have explicitly 
checked with the PYTHIA Monte Carlo results that our 
interpolation is accurate at the percent level both on 
the reconstructed pion yield and on the final gamma-ray 
spectrum from tt" decay, as given by Eq. 

The contribution to the 7-ray spectrum from produc- 
tion and decay of mesons other than pions (mostly 77, 
77', charmed and bottom mesons) and of baryons is usu- 
ally subdominant as compared to tt'^ decay. These ad- 
ditional contributions can be safely neglected (they typi- 
cally contribute only up to 10% of the total flux for < 
1 GeV). A notable exception is given by the hadroniza- 
tion of bb pairs at low production energies, i.e. for neu- 
tralino masses between the production threshold for a 
&"meson and about 10 GeV. In this case jet flavor con- 
servation leads to the production of a bottom meson 
B = S", 5° with 100 % probability. In the PYTHIA 
code, 75% of the times the B meson is in the excited 
state and decays through B* B + "/ with mg — ms — 
46 MeV. Since the B* mesons are produced almost at 
rest (to_b — 5.3 GeV) for < 10 GeV, they generate 
a (slightly boosted) gamma-ray line that dominates the 
other contributions below E^ ~ 100 MeV. We have then 
included this peculiar contribution to our interpolating 
procedure fjli] . 

Neutralino annihilation into Icpton pairs can also pro- 
duce gamma-rays from electromagnetic showering of the 
final state leptons. This process can be dominant for 
Ej < 100 MeV, when the neutralino annihilation pro- 
cess has a sizable branching ratio into lepton pairs. In the 
case of production of r leptons, their semihadronic decays 
also produce neutral pions, which then further contribute 
to the gamma-ray flux. Also these additional contribu- 
tions arc included in our numerical evaluations, again by 
modelling the gamma-ray production with the PYTHIA 
Monte Carlo for the same set of neutralino masses quoted 
above, and by numerically interpolating for other values 
of m^. 

When the neutralino masses are sufficiently large, the 
annihilation channels into Higgs bosons, gauge bosons 
and tt pairs become kincmatically accessible. We com- 
pute analytically the full decay chain down to the pro- 
duction of a quark, gluon or a Icpton. The ensuing 7-ray 
spectrum is then obtained by using the results discussed 
above for quarks, gluons and leptons, by properly boost- 



ing the differential energy distribution to the rest frame 
of the annihilating ncutralinos (see e.g. Appendix I in 
Ref. ^ for details). 

B. The geometrical factor 

The integral along the line of sight in Eq. © is 
the quantity that takes into account the shape of the dark 
matter profile. For small values of "0, /(■0) is very sen- 
sitive to possible enhancement of the density at the GC 
and therefore large differences for the different density 
profiles are expected. 

When comparing to experimental data, Eq. jSJ is aver- 
aged over the telescope aperture-angle Atp: 

Im = ^I mdi>. (7) 

The gamma-ray fiux is therefore proportional to /ai/'- 

In our analysis on the galactic center emission, we will 
employ data from the EGRET experiment |33ll3q . whose 
angular resolution is given by the longitude-latitude 
aperture: |A/| < 5°, |A6| < 2°. Table IJ shows the val- 
ues of /ai/j for the density profiles discussed in Sec. IIIII 
The effect of changing the core radius of an isothermal 
sphere is not negligible: an increase of a factor of 2.3 is 
obtained by reducing a from 3.5 kpc to 2.5 kpc. In the 
case of singular distributions a small-distance cut-off of 
Tc = 0.01 pc is assumed (inside Tc the density is assumed 
to be constant). A NEW profile then gives a fiux which 
is about 10 times larger than an isothermal sphere with 
a = 3.5 kpc. The very steep Moore et al. profile would 
produce a fiux about 60 times larger than a NEW. The 
Navarro et al. profile with r-dependent log-slope sits 
between the NEW and Moore et al. cases: though not 
singular, it nevertheless provides a fiux which is about 
3 times larger than a singular NEW halo. The variabil- 
ity of /a^, on the dark matter profile can therefore be as 
large as a factor of 600, comparing the very steep Moore 
et al. profile with an isothermal sphere with a large core 
radius. However, the recent critical analysis on numeri- 
cal simulation of Ref. implies that a factor of 30 is 
likely to be a more plausible interval. 

In the case of high galactic latitudes, the dependence of 
/ai/j on the density profile is much milder. At these high 
latitudes, EGRET identifies a residual gamma-ray flux 
which is ascribed to a possible extragalactic component, 
but could as well be due to dark matter annihilation in 
the galactic halo. We will consider in our analysis two 
different angular regions: |6| > 10° with the exclusion 
of |;| < 40° and 10° < |5| < 30° around the galactic 
center "fi^ (region A) and \b\ > 86° (region B). Region 
B has been considered in the rcanalysis of EGRET data 
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Isothermal 


Isothermal 


NFW 


Moore et al. 


r-dependent log-slope Eq.JJJ 


a = 3.5 kpc 


a = 2.5 kpc 


a = 25 kpc 
Tc = 0.01 pc 


a = 30 kpc 
rc = 0.01 pc 


Q = 0.142 
r_2 = 26.4 kpc 
P-2 = 0.035 GeV cm"^ 


18.5 


42.5 


184.2 


10866 


600 



TABLE I: Values for Ja^ in Eq.Q when different dark matter distributions are assumed (in units of GeV^ cm~^ kpc). The 
angular region of integration is defined by the intervals; |AZ| < 5°, |A&| < 2°. The first two columns refer to an isothermal 
distribution with a core a = 3.5 and 2.5 kpc, respectively. The third and fourth columns refer to singular DM distributions: a 
NFW with a scale length a = 25 kpc and a Moore et al. profile with a scale length a = 30 kpc; in both cases, the DM profile 
has a cut-off radius rc = 0.01 pc. The last column refers to the density profile of Eq.Q with the parameters of the distribution 
Gl in Table III of Ref. [H: a = 0.142, r_2 = 26.4 kpc and p-2 = 0.035 GeV cm~^ For all these profiles: pi = 0.3 GeV cm'^ 



by the authors of Ref. j3a | , where more stringent limits 
on the the gamma-ray residual intensity from the galac- 
tic poles have been derived. The value of /ai/j for region 
A is 1.66 for a NFW profile, and ranges from 1.61 for 
the isothermal sphere with a = 3.5 kpc to 1.80 for the 
r-dependent log-slope profile: when pointing away from 
the critical behavior of the density profiles at galactic 
center, the line of sight integral is almost universal. In 
the case of region B, /a^ = 0.67 for the NFW distri- 
bution, and ranges from 0.62 for the isothermal sphere 
with a = 3.5 kpc to 0.69 for the r-dependent log-slope 
profile. In both cases, the Moore et al. profile gives a 
line-of-sight integral slightly smaller than in the case of 
the r-dependent log-slope profile. 

C. Signal from the galactic center 

Data from low galactic latitudes < 10°), includ- 
ing the galactic center region, have been collected by the 
EGRET telescope The diffuse gamma-ray flux of 

the inner galaxy measured by EGRET shows a possible 
excess over the estimated background at energies larger 
than about 1 GeV. 

Clearly a firm assessment of an excess requires a good 
knowledge of the standard production of gamma rays in 
our Galaxy. At the energies of interest for our analy- 
sis - namely, from about 100 MeV to tens of GeV - the 
main production mechanism of 7-rays is the interaction 
of cosmic rays (mainly protons and helium nuclei) with 
the interstellar medium (atomic and molecular hydrogen, 
and helium). In these strong reactions tt^'s are produced, 
and hence 7-rays via pion decay: tt'^ 27. The ensu- 
ing spectrum has a bump around 70 MeV, and drops at 
high energies with an energy power law which follows the 
progenitor cosmic ray spectrum {E~", with a ~ 2.7). 
Another source of 7-rays comes from inverse- Compton 
scattering of cosmic ray electrons off the interstellar pho- 
tons. In particular, energetic electrons may scatter off 
the cosmic microwave background, and off the infrared. 



optical and ultraviolet radiation arising from stellar ac- 
tivity and dust. The third radiation component origi- 
nates from electron bremsstrahlung over the interstellar 
medium, which may be partially or even totally ionized. 
Bremsstrahlung 7-rays are mostly important in the low 
energy tail. For a full calculation of these three main ra- 
diation components one needs a good knowledge of the 
physics of cosmic rays and of the interstellar medium in 
the region of interest. This is particularly unlikely when 
dealing with the galactic center area. 

In the literature, several different results have been 
achieved on the subject. First of all, the EGRET Col- 
laboration developed a detailed calculation of the 7-ray 
background at the energies of interest for the detector 
[33L W\ : this calculation shows a clear deficit of 7-rays 
towards the GC with respect to measurements. The ex- 
cess in the data is apparent for > 1 GeV, where the 
shapes of the spectra of the estimated background and 
the data differ significantly. At lower energies the spec- 
tral agreement is instead rather good. Similar conclu- 
sions have been drawn in Refs. |33| with some different 
procedure in the calculation of the background. In this 
paper a harder, probably unrealistic, nucleonic spectrum 
is shown not to be anyway sufficient to explain the GeV 
excess. Some modifications towards harder electron and 
nucleon spectra are studied in Ref. |39| . but a satisfac- 
tory agreement with data is not achieved (notice that it 
is anyway hard to reconcile these hypothesis with galac- 
tic cosmic rays measurements). The results of all these 
analyzes favor the interpretation of the EGRET data in 
terms of an excess over the background, mostly at ener- 
gies above 1 GeV. 

However, different assumptions on acceleration |4o| 
and diffusion |4 j | of cosmic ray nucleons, and on the spec- 
tral shape of primary nucleons in the interstellar space 
[4^ , have been proposed and lead to a quite good agree- 
ment of the EGRET measured flux: almost all the spec- 
tral features are reproduced by these calculations. In this 
case, the EGRET data would be explained in terms of 
the standard galactic 7-ray production. 
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FIG. 1: Gamma-ray flux from the galactic center inside the an§u?ar region \Al\ < 5°, |A6| < 2° for a NFW matter density 
profile. The scatter plots are derived by a full scan of the parameter space of the supersymmetric model described in Sec. HTl 
Crosses (red) and dots (blue) denote neutralino configurations with 0.095 < Q-^h^ < 0.131 and Q,-^h^ < 0.095, respectively. Top 
LEFT: flux calculated at Ej = 0.12 GeV; the horizontal line shows the gamma-ray flux measured by EGRET 1331 . assumed to 
be compatible with the estimate of the background |33|. Top right: flux calculated at Ej = 1.5 GeV; the solid horizontal 
line shows the gamma-ray flux measured by EGRET l33ll . the dashed line is an estimate of the gamma-ray background l33l . 
Bottom: flux calculated at E^y = 15 GeV; the solid horizontal line shows the gamma-ray flux measured by EGRET [j^l . the 
shaded horizontal band denotes the la error bar on the EGRET data and the dashed line is an estimate of the gamma-ray 
background [s^. 



At present, it is very difficult to favor one model 
against the others, both on theoretical and observational 
basis. This implies that the uncertainty on the calcula- 
tion of the galactic 7-ray flux, and in particular at the 
galactic center, is very difficult to quantify. 

Due to these open problems on the determination of 
the background component, we will develop our analysis 



along two paths. First of all we will discuss whether, and 
under which conditions, it is possible to set constraints on 
low-mass neutralinos from the gamma-ray studies. Then 
we will comment on the possibility for low-mass neutrali- 
nos to explain the putative EGRET excess. 

The gamma-ray flux from the galactic center inside the 
angular region |A;| < 5°, |A6| < 2° for a NFW matter 
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FIG. 2: Gamma-ray spectra $^(_E^), multiplied by E?f, from the galactic center inside the angular region \Al\ < 5°, |A6| < 2°, 
as functions of the photon energy. Left panel: the dotted line is the spectrum for a neutralino with mass = 30 GeV, 
calculated for a density profile with a factor 30 of enhancement with respect to the NFW case; the dashed line is the gamma ray 
background calculated in Ref. [j^l . reduced by 10%; the solid line is the total flux, sum of the supersymmetric signal and the 
background; the experimental points are the EGRET data [s^- Right panel: the same, for = 40 GeV and for a density 
profile with a factor 32 of enhancement with respect to the NFW case. Both supersymmetric configurations have been selected 
from the points shown in Fig. The numbers quoted in the legend inside parentheses denote the values of the neutralino 
annihilation branching ratios into bb and fr. 



density profile is shown in Fig. ^ at three representa- 
tive photon energies: = 0.12, 1.5 and 15 GeV. These 
energies correspond to three energy bins of the EGRET 
detector. The scatter plots of the two top panels display a 
peculiar funnel shape at small masses. This is due to the 
fact that the neutralino flux is bounded from below by the 
cosmological limit fl^h^ < (r2c-DM/i^)max- This feature 
is similar to the one wc found in Refs. [ij, |3j , in con- 
nection with neutralino direct detection rates. The vari- 
ation in shape of the scatter plots, when is increased, 
is easily understood in terms of Eq.0). At E^ ~ 0.12 
GeV the m~'^ behavior is clearly visible. Energies of the 
order of 100 MeV are very crucial in offering the pos- 
sibility to set limits on the very light neutralino sector. 
By increasing the photon energy, the lightest ncutrali- 
nos do not have enough phase space to produce photons 
at this energy (since they annihilate almost at rest in 
the Galaxy): therefore the gamma-ray flux at very low 
masses becomes progressively depressed, as E^ increases. 
At E^ — 1.5 GeV a neutralino with a mass of 6 GeV can 
produce approximately the same flux as a 10-15 GeV neu- 
tralino. At Ey ~ 15 GeV all the neutralinos lighter than 
15 GeV obviously do not produce any photon: at this 
energy the maximal fluxes are obtained for neutralinos 
with masses around 30 GeV. 

The first conclusion which can be drawn from Fig. ^ 
is that the supersymmetric model considered here is not 



constrained, at present, by EGRET data for a NFW den- 
sity profile. An increase in the flux by a factor of 3.3, as 
would be in the case of the 7'-depcndent log-slope profile 
of Eq.©, is also not enough to set limits. 

Larger enhancements in the geometrical factor /a^ 
over the NFW case are necessary to set limits. The com- 
parison of the three panels of Fig. ^shows that the limits 
come from different energies E^, depending on the neu- 
tralino mass. For very light neutralinos, i.e. < 10 
GeV, the lowest energy bin is the relevant one. In this 
case an enhancement of a factor of 6 would allow to raise 
the predicted fluxes for E~^ ~ 0.12 GeV at the level of 
the EGRET measurement, and therefore to start setting 
limits. For masses in the range 10 GeV < niy^ < 20 GeV 
the E~^ = 1.5 GeV bin sets more stringent limits, at least 
on a fraction of the supersymmetric models, but only for 
a factor of enhancement of at least 15-20 over the NFW 
case. These factors are pretty large, even though not 
as large as the one which refers to a Moore et al. pro- 
flle, which is about 60, as discussed before. For masses 
around 30-40 GeV the best limits come from the highest 
energy bin i?-, = 15 GeV, where a factor of 20-25 would 
allow the fluxes to reach the EGRET data. In the case 
of the standard MSSM, where the neutralino has masses 
larger than 50 GeV the lowest energy bins are always less 
constraining that the higher energy ones, as can be seen 
by comparing the different panels of Fig. Instead, the 



9 



lower energy bins are crucial for the study of the low- 
mass neutralinos. We finally comment that a Moore et 
al. profile would make all the fluxes for < 10 GeV 
incompatible with the data, but this profile is less likely, 
as we discussed above. 

Now, let us turn to a brief discussion of the possibil- 
ity to explain the EGRET excess by means of low-mass 
neutralinos. The analysis made above on the behavior 
of the gamma-ray fluxes in the three representative en- 
ergy bins of Fig. ^ shows that neutralinos in the mass 
range 25 GcV < < 40 GcV are the ones which may 
have the possibility to fill the excess in the energy range 
above 1 GeV, without spoiling the lower energy behavior 
of the background which is supposed to have an accept- 
able agreement with the data. We show that indeed the 
low-mass neutralinos in this mass range are able to ex- 
plain the EGRET excess in Fig. |2 In this figure we 
plot the predicted gamma-ray spectra for two represen- 
tative supersymmetric configuration when, for definite- 
ness, the gamma-ray background as calculated in Ref. 
[srt is assumed. In the left panel of Fig. [3 we show a 
supersymmetric configuration with m-^ ~ 30 GcV and a 
relic abundance in the proper range to explain dark mat- 
ter: fl^h^ = 0.12. The dominant annihilation channels 
of these low-mass neutralinos are xx ~^ t't &iid XX ^ bb 
P, Gamma rays coming from annihilation into r's 
give a harder spectrum as compared to the b channel. 
In this representative point the two channels have (ap- 
proximately) the same branching ratio: this gives siz- 
able contribution to the gamma-ray fiux in all the en- 
ergy range from 1 to 10 GeV, which is where the excess 
in the EGRET data is more pronounced. By allowing 
the background to fluctuate down by 10% and by using a 
geometrical factor /ai/j 30 times larger than in the NEW 
case, we show that a pretty good agreement between the 
total flux and the data can be obtained. We are not 
quoting a statistical significance for this agreement since 
we are not performing a systematic statistical analysis 
here: however, we are interested in showing that, in ad- 
dition to the standard MSSM neutralinos with masses 
larger than 50 GeV 0] , also low-mass neutralinos have 
the capability of explaining the putative EGRET excess. 
In both cases, low-mass and standard neutralinos, the 
values of the line-of-sight integrals Ia^p which are able 
to explain the EGRET excess are much larger than what 
is provided by a NFW density profile. However, for neu- 
tralinos in the 30-40 GeV mass range these enhancement 
factors are smaller than in the case of heavier neutralinos, 
due to the behavior. 

The right panel of Fig. |21 shows a second example, 
with a neutralino of = 40 GeV and fJ^ft,^ = 0.11. 
The background of Ref. [s^l is again scaled down by 
10%. In this case the geometrical factor is 32 times the 



NFW one. The branching fraction of annihilation into b 
quarks is larger than in the previous case. This enhances 
the contribution to the gamma-ray flux in the 1-3 GeV 
range without spoiling the agreement at larger energies. 

A detailed analysis of the spectral features of the 
gamma-ray fluxes produced by neutralino annihilation 
and their comparison with the EGRET data is beyond 
the scope of the present paper and will be presented else- 
where. 



D. Signal from high galactic latitude 

Data from high galactic latitudes have been collected 
by the EGRET telescope . An analysis of the mea- 
surements taken over the latitudes |6| > 10°, and exclud- 
ing the region |/| < 40° and 10° < |6| < 30° around the 
galactic center, has been performed in Ref. [i^l- AH the 
identifled sources as well as the components due to the in- 
teractions of cosmic rays with the galactic disk gas have 
been subtracted [sJl- The residual flux, averaged over 
the considered portion of the sky, has been shown to be 
isotropic and well fltted by the power law ^fjj^^'^ (Ej) = 
kiE^/Eo)-", where k ^ (7.32±0.34)-10-s photons cm-^ 
scc-i sr-1 GeV-i, a = 2.10 ± 0.03 and Eq = 451 MeV. 
This spectrum is often referred to as the extragalactic dif- 
fuse emission, since no known source inside the Galaxy 
seems to be responsible for it. One possibility is that it 
is due to unresolved gamma-ray-emitting blazars. Rely- 
ing on the analysis by EGRET, one can use the residual 
flux as an upper bound to any flux due to exotic sources, 
including annihilation of relic neutralinos. Recently, a 
re-analysis of EGRET data has been performed in Ref. 
[iif . taking particular care of the spatial dependence of 
the observed photons. Working on the integrated flux 
and taking into account contributions from several galac- 
tic tracers, the authors of Ref. [s^ show that the high- 
latitude 7-ray sky exhibits strong galactic features and 
could be well accommodated by simple galactic models. 
Conservative constraints have been set on the flux inte- 
grated above 100 MeV and averaged over different sectors 
of the sky far from the galactic plane. In this scenario, 
the room left to an unexplained diffuse flux - often con- 
sidered as an extragalactic background, but also possibly 
due to exotic galactic sources - is much smaller (by a 
factor of three, at least) than the one reported in Ref. 
[i^ l. Here we consider the upper limit of Ref. [s^ on 
this possible residual, isotropic flux in the polar re- 
gion (|6| > 86°): < 0.6 ■ 10"^ 7 sec-^cm-^sr-i, and 
compare it to our estimates for the 7-ray flux due to neu- 
tralino annihilation averaged in the same spatial region. 

The results for both estimates are shown in Fig. Ofor 
a NFW profile. As in the case of the galactic center emis- 
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< 0.095, respectively. Left panel: fiux calculated at = 0.12 GeV in 



the high-latitude regions defined by \b\ > 10° with the exclusion of |/| < 40° and 10° < |&| < 30° around the galactic center; 
the horizontal line shows the gamma-ray residual flux identified by EGRET 0|. Right panel: integrated fiux for energies 
above = 0.1 GeV in the polar regions defined by |b| > 86°; the dashed horizontal line shows the upper limit on a possible 
residual flux in the polar regions, obtained in Ref. |35l . 



sion, also for high latitudes we do not have constraints 
on the supersymmetric parameter space. Contrary to the 
case of the galactic center region, for high latitudes the 
geometrical factor /a^ is practically independent of the 
halo profile, as discussed before. We therefore conclude 
that at present the 7-ray signals from high galactic lati- 
tudes do not provide any constraint on the supersymmet- 
ric parameter space. The situation could change if fur- 
ther studies will show that a much bigger fraction of the 
EGRET measured flux at high latitudes is due to galac- 
tic foreground or when next-generation experiments will 
provide further information. 

V. ANTIPROTONS IN COSMIC RAYS 

Like in the case of the gamma-ray flux, the produc- 
tion of antiprotons from neutralino annihilation results 
from the hadronization of quarks a nd g luons created in 
the annihilation process [sa, H [43 . The differ- 

ential rate per single annihilation, unit volume and time 
is defined as: 

,r'-ir.T,)^i.^.,..)'^l[^f. (8) 

where Tp denotes the antiproton kinetic energy, dNp /dTp 
(indicated as g{Tp) in Refs.js^^^) is the differential an- 
tiproton spectrum per annihilation event, and the factor 



1/2 accounts for the number of annihilating neutralino 
pairs. The spectrum dNp/ dTp is evaluated by means 
of the Monte Carlo simulations we already used in Sec. 

evil 

Once antiprotons are produced in the dark halo, they 
diffuse and propagate throughout the Galaxy. To de- 
scribe these processes, we follow the treatment of ReL 
[3^ . to which we refer for details. Here we only recall that 
the propagation of antiprotons has been considered in a 
two- zone diffusion model 0, 0, [s^l , defined in terms of 
six free parameters whose role is to take into account the 
main physical processes which affect the propagation of 
cosmic rays in the Galaxy: acceleration of primary nuclei, 
diffusion, convective wind, reacceleration process and in- 
teraction with the interstellar medium. These free pa- 
rameters are constrained by comparing the fluxes of var- 
ious cosmic ray species calculated in our diffusion model 
with observations. In this regard, the most important 
observable is the measured Boron/Carbon ratio {B/C), 
whose analysis within our diffusion model is presented in 
Ref. ji^. The parameters constrained by the B/C mea- 
surements have been shown to be compatible with a series 
of other observed species 0, |^ [H^l , further supporting 
the employed model. Therefore, in the calculation of the 
primary antiproton flux we only use those values for the 
propagation parameters which provide a good statistical 
agreement to the B /C data. One of the main results ob- 
tained in Ref. |32l | is that the supersymmetric antiproton 
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FIG. 4: Antiproton flux at Tp = 0(23 GleV as a function of the neutralino mass, calculated at solS,r minimum. The scatter plots 
are derived by a fuU scan of the parameter space of the supersymmetric model described in Sec. EI ^ spherical isothermal 
dark matter density profile has been used. The solar modulation is calculated at the phase of solar minimum. Crosses (red) 
and dots (blue) denote neutralino configurations with 0.095 < ^^h?' < 0.131 and Q^h^ < 0.095, respectively. The shaded 
region denotes the amount of primary antiprotons which can be accommodated at Tb = 0.23 GeV without entering in confiict 
with the experimental BESS data Im. Is^ and secondary antiproton calculations [sif . Left panel: the best fit set for the 
astrophysical parameters is used. Right panel: the astrophysical parameters which provide the most conservative antiproton 
fluxes are used. 



flux, when calculated with the selected propagation pa- 
rameters, is affected by a large uncertainty. At low energy 
this uncertainty reaches almost two orders of magnitude, 
while it diminishes to a factor of thirty at higher ener- 
gies. Only better and more complete data on cosmic rays 
(both stable and radioactive) could help in reducing this 
uncertainty. 

At variance with the gamma ray signal, the supersym- 
metric antiproton flux measured at Earth is almost in- 
sensitive to the specific form of the dark matter distribu- 
tion function, among those discussed in Sec. IIIII Indeed, 
these various distributions differ mainly at the galactic 
center, while at solar neighborhood differences are very 
mild. Since charged particles, such as antiprotons, suf- 
fer enormous energy redistributions, gains and losses, it 
is almost impossible for an antiproton produced around 
the galactic center to reach the Earth. This property 
was shown in Ref. |53| . and quantified in Ref. |32j for 
the case of a NEW distribution and an isothermal one. 

In the present work, we use directly the results reached 
in Ref. 1321 . where the function 



supersymmetric flux factor: 



T5(Tp) 



(9) 



was calculated. In this equation, ^p{rQ,Tp) is the in- 
terstellar antiproton flux after propagation and T is the 



T 



1 



(10) 



The quantity C^nsyC^p) may be considered as a measure 
of how the source flux q^^^{r, Tp) is modified by the prop- 
agation of antiprotons in the Galaxy before reaching the 
heliosphere. In the results presented in the following, we 
have calculated the antiproton flux <^p(rQ, Tp) according 
to Eq. ©, where the C^™^(Tp) function has been taken 
directly from Ref. for a few representative combina- 
tions of the propagation parameters and source spectra 
g{Tp). We have calculated the quantities entering the 
factor T as described in Sec. |ll| 



A. Secondary antiprotons 

Antiprotons in the Galaxy are also produced via stan- 
dard interactions. Proton and helium cosmic rays in- 
teract with the interstellar hydrogen and helium nu- 
clei, producing quarks and gluons that subsequently can 
hadronize into antiprotons. A calculation of this sec- 
ondary antiproton flux has been done in Ref. |5lj . to 
which we refer for details. Here we only emphasize the 
main results obtained in that work: i) The antiproton 
flux has been evaluated consistently by employing the 
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propagation parameters as derived from a full and sys- 
tematic analysis on stable nuclei 0| ; ii) This secondary 
antiproton flux is in very good agreement with the data 
taken from the experiments BESS Is^. [s^, AMS [s^, 
CAPRICE (see Fig. 14 in Ref. [sj); iii) The uncer- 
tainty on the final flux due to propagation is about 20%, 
with a slight dependence on the energy. Another impor- 
tant source of uncertainty of order 20%-25% resides in 
the nuclear production cross sections, in particular when 
considering the interactions over the interstellar helium. 



B. Constraints on a primary antiproton source? 

As discussed above, the secondary antiproton flux al- 
ready provides a satisfactory agreement with current ex- 
perimental data, and then no much room is left to pri- 
mary contributions. This situation suggests that antipro- 
ton data could be used to place significant constraints on 
supersymmetric parameters. However, one has to notice 
that, as shown in |32l |. the supersymmetric primary flux is 
affected by uncertainties much larger than those related 
to the secondary flux. This is due to the fact that the 
sources of the latter are located in the galactic disk. On 
the contrary, the relic ncutralinos are expected to be dis- 
tributed in the whole galactic halo and then produce an 
antiproton flux much more sensitive to the astrophysical 
parameters. 

To show quantitatively how the experimental data 
could constrain the supersymmetric parameters, in Fig. 
01 we display the antiproton flux evaluated at Tp = 0.23 
GeV for a full scan of our supersymmetric model de- 
scribed in Sec. EI As expected, the scatter plot is promi- 
nent at small masses. Furthermore, it is remarkable that 
for nix ^ 25 GeV the scatter plot is funnel-shaped. The 
reason is the same as the one given above in connection 
with Fig. n The two panels of Fig. 0] correspond to two 
different sets of the propagation parameters. One, used 
in the left panel, is the set giving the best fit to the B/C 
ratio, while the other, hereby denoted as the conservative 
set and used in the right panel provides the lowest (sec- 
ondary and primary) antiproton fluxes. A spherical cored 
isothermal distribution for dark matter has been used. 
However, as mentioned before, a different choice does 
not significantly modifies the scatter plots. The shaded 
region denotes the amount of primary antiprotons which 
can be accommodated at Tp = 0.23 GeV without enter- 
ing in conflict with the BESS experimental data [s^lssl 
and secondary antiproton calculations |5lj . 

From the right panel of Fig. ^wc conclude that, within 
the current astrophysical uncertainties, one cannot de- 
rive any constraint on the supersymmetric parameters, 
if one assumes a very conservative attitude in the selec- 



tion of the propagation parameters. It is worth noticing 
that even within this choice, some supersymmetric con- 
figurations at very small masses are close to the level of 
detectability. As a further comment on the left panel of 
Fig. 21 we wish to stress that any further breakthrough 
in the knowledge of the astrophysical parameters would 
allow a significant exploration of small mass configura- 
tions, in case the conservative set of parameters is ex- 
cluded. Should the effect of antiproton propagation turn 
out to be equivalent to the one obtained with the best fit 
set, the analysis of cosmic antiprotons would prove quite 
important for exploring very light neutralinos. This is 
particularly true for neutralino masses below < 15 GeV, 
in view of the typical funnel shape displayed in the scat- 
ter plots. 

VI. UPGOING MUONS AT NEUTRINO 
TELESCOPES 

Indirect evidence for WIMPs in our halo may be ob- 
tained at neutrino telescopes by measurements of the 
upgoing muons, which would be generated by neutri- 
nos produced by pair annihilation of ncutralinos captured 
and accumulated inside the Earth and the Sun |58l l59l| . 
The process goes through the following steps: capture 
by the celestial body of the relic neutralinos through a 
slow-down process due essentially to neutralino elastic 
scattering off the nuclei of the macroscopic body; accu- 
mulation of the captured ncutralinos in the central part 
of the celestial body; ncutralino-ncutralino annihilation 
with emission of neutrinos; propagation of neutrinos (we 
have included the v^j.-Vr vacuum oscillation effect with 
parameters: Ato^ = 3 x 10""^ eV^, sin 6' = 1) and con- 
version of their i'^ component into muons in the rock 
surrounding the detector; propagation and detection of 
the ensuing upgoing muons in the detector. 

The various quantities relevant for the previous steps 
are calculated here according to the method described 
in Rcfs. [i^, to which we refer for further details. We 
just recall that the neutrino flux due to the annihilation 
processes taking place in a distant source like the Sun, 
as a function of the neutrino energy E^, is given by 

dE, And^^^ dE, ' ^ ' 

where Ta is the annihilation rate inside the macroscopic 
body [60|, d is the distance from the source, F de- 
notes the x~X annihilation final states, B^j denotes the 
branching ratios into heavy quarks, r leptons and glu- 
ons in the channel F; dNf^/dE^ is the differential dis- 
tribution of the neutrinos generated by the hadroniza- 
tion of quarks and gluons and the subsequent hadronic 




FIG. 5: Flux of upgoing muons aS a funcAon of the neutralino mass. The scatter plots are derived oy a, scan of the parameter 
space of the supersymmetric model described in Sec. [H] Crosses (red) and dots (blue) denote neutralino configurations with 
0.095 < fl^h^ < 0.131 and Q^h^ < 0.095, respectively. Left panel: signal from the Earth; the solid, dashed and dotted lines 
denote the experimental upper limits from SuperKamiokande [ij, MACRO [i^ and AMANDA [g^, respectively. Right panel: 
signal from the Sun; the solid, dashed, dot-dashed and dotted lines denote the experimental upper limits from SuperKamiokande 
[e^ . MACRO lil, Baksan [tI and AMANDA respectively. 



semileptonic decays. The annihilation rate is given by 
Fa = C/2 tanh^ {to / ta) where to is the age of the 
macroscopic body, ta = {CCa)~^^'^ , Ca is the annihila- 
tion rate proportional to the neutralino-neutralino anni- 
hilation cross-section and C denotes the capture rate per 
effective volume of the body. The capture rate C is calcu- 
lated here as in Refs. [illliil, the WIMP velocity distri- 
bution in the galactic inertial frame being described by a 
Maxwellian function with a dispersion velocity of 270 km 
s~^. We recall that the capture rate by the Earth is fa- 
vored when the WIMP mass is close to the nuclear mass 
of one of the main chemical components of the Earth: 
Oxygen, Magnesium, Silicon in the mantle and Iron in 
the core [g^I . We have neglected the contributions of the 
light quarks directly produced in the annihilation process 
or in the hadronization of heavy quarks and gluons, be- 
cause these light particles stop inside the medium (Sun 
or Earth) before their decay. For the case of the Sun we 
have also considered the energy loss of the heavy hadrons 
in the solar medium. 

A. Neutrinos from the center of the Earth 

In the left panel of Fig. [Slwc show the expected upgo- 
ing muon flux integrated for E^^ >1 GeV, as a function 
of m^, and compared to the present experimental upper 
bounds on the same quantity from the experiments Su- 
perkamiokande MACRO 0, and AMANDA 



For < 40 GcV the signal from the Earth presents sev- 
eral peaks due to resonant capture on Oxygen, Silicon 
and Magnesium (we recall that the process of capture 
on Earth is driven by the coherent neutralino-nuclcus 
cross section). These elements are almost as abundant 
in Earth as Iron, which is the most important target nu- 
cleus for neutralino capture at higher masses. The dip at 
~ Mz/'i- is due to the rescaling prescription of Eq. 
(13, since the resonance in the Z-exchange annihilation 
cross section reduces the relic abundance Q,^h? . More- 
over, for < 25 GeV, the branching ratio B^J to the 
rf final state, which is the one with the highest neutrino 
yield per annihilation, is suppressed. This last property 
is due to the fact that, in this range of m^, the final state 
to bb in the annihilation cross section is required to be the 
dominant one in order to keep the relic abundance fl^^h^ 
below its cosmological upper bound Q. This, together 
with the fact that lower x masses imply softer v spec- 
tra that produce less ^'s above threshold, explains why 
the up-going muon signal expected for light neutralinos 
(m^ < 50 GeV) turns out to be below the level reached 
at higher masses. 

It is worth noting that substantial modifications to 
the standard Maxwellian velocity distribution in the so- 
lar neighborhood have been considered in recent years, 
with two confiicting models. Damour and Krauss 
proposed the existence of a solar-bound population, gen- 
erated by WIMPs which scattered off the Sun surface 
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and were then set, by perturbations from other planets, 
into orbits crossing the Earth, but not the Sun. The 
existence of this low-speed population would make the 
capture of WIMPs by the Earth very efficient with an en- 
suing dramatic enhancement of the up-going muon flux 
as compared to the standard case . On the contrary, 
Gould and Alam used arguments based on calcula- 
tions of asteroid trajectories to conclude that solar-bound 
WIMPs could evolve quite differently, inducing a signif- 
icant suppression in the up-going muon flux from the 
Earth, as compared to the standard Maxwellian case, for 
WIMP masses larger than 150 GeV. A very recent re- 
analysis of this problem |69| supports the conclusions of 
Gould and Alam, though with a less dramatic suppres- 
sion effect. Our results have been derived using the stan- 
dard Maxwellian distribution. The results of Ref. 
would not significantly alter our conclusion on the de- 
tcctability of light ncutralinos. The maximal fluxes are 
obtained for resonant capture on O, Si and Mg nuclei in 
the mantle: in this situation no suppression occurs. For 
ncutralino masses away from the resonant condition, a 
reduction factor up to 0.8 for mo < < mpe and up 
to 2-3 for < 10 GeV is possible : however, in these 
cases, the upgoing muon flux is already very depressed, 
as it is shown in Fig. El left panel. In conclusion, using 
the standard Maxwellian distribution, the present mea- 
surements of up-going muons from the Earth put some 
constraints on ncutralino conflgurations for masses above 
50 GeV. For lighter neutralinos, explorations by neutrino 
telescopes would require a substantial increase in sensi- 
tivity while keeping a low energy threshold (close to 1 
GeV). This in turn would imply a sizable extension of the 
telescope and a dense array of photomultipliers, which is 
certainly feasible, but very expensive. 



B. Neutrinos from the Sun 

In the right panel of Fig. [S] we show the up- 
going muon flux expected from the Sun, integrated for 
>1 GeV, as a function of m^. The signal is com- 
pared to the present experimental upper bounds on this 
flux coming from the experiments Superkamiokande|63l| . 
MACRO [ej, AMANDAjii and Baksanlz^. Also in 
this case the signal level turns out to be suppressed for 
™x < 50 GeV as compared to what is obtained at higher 
masses, the reasons of this behavior being the same as in 
the case of the Earth. On the other hand, the enhance- 
ment of the signal at mw is due to a peculiar 
behavior of the neutralino-nucleon spin-dependent cross 
section, which drives the ncutralino capture in the Sun 
(mainly on hydrogen). This cross section reaches its max- 
imum whenever the Z-cxchangc channel dominates, and 



this requirement is verified when the neutralino-Z cou- 
pling, proportional to the combination a\—a\, is maximal 
|7ll |. By numerical inspection we have verified that this 
last quantity is significantly peaked for ~ M^/. In 
this range of masses the annihilation channel to W~^W~ 
opens up and dominates the annihilation cross section. 
We conclude here that investigations of light neutrali- 
nos by up-going muons from the Sun do not provide 
favourable prospects. 



VII. CONCLUSIONS AND PERSPECTIVES 

We have considered the most relevant indirect strate- 
gies for detecting the presence of relic neutralinos in our 
Galaxy through the products of their self-annihilation. 
This includes annihilations taking place directly in the 
galactic halo or inside celestial bodies (Earth and Sun). 
Our investigation has been performed in the frame of an 
effective supersymmetric model at the electroweak scale 
with no assumption on gaugino mass universality at the 
GUT scale. The range of the ncutralino mass taken into 
consideration brackets a wide interval, from 6 GeV up to 
500 GeV. While the low extreme is decided by the lower 
bound of 6 GeV, estabhshed in Ref. 0,13; the upper ex- 
treme of 500 GeV is chosen only for convenience. Actu- 
ally, no model-independent upper limit for the ncutralino 
mass is available, apart for a generic value of order of 1 
TeV, beyond which the raison d'etre of supersymmetry 
fades away. Though our calculations span over the wide 
range of the ncutralino masses recalled above, our discus- 
sions were focused on light neutralinos, i.e. neutralinos 
with masses below 50 GeV: this value corresponds to the 
lower bound of when gaugino-mass unification is as- 
sumed. Indeed, light neutralinos are rarely considered 
in the literature, though their properties are quite in- 
teresting, as we already proved in connection with their 
cosmological properties and their detectability by cur- 
rent experiments of direct WIMP search Q, 01 ■ Thus 
the present paper is the natural continuation of our previ- 
ous investigations on neutralinos of small mass. Different 
galactic dark matter distributions have been considered, 
from the cored isothermal one to the profiles obtained 
by numerical cosmological simulations in Refs. |l6l Il7l |. 
including also the most recent ones of Ref. ^3 ■ 

Now we summarize the main results of the present pa- 
per: 

• For the 7 rays we have considered separately fiuxes 
from the galactic center and high latitude regions, 
and compared our predictions with the EGRET 
data. Our numerical results have been provided 
employing as a reference DM distribution the NFW 
profile. We have shown that in this case the 
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EGRET data at all angles do not put any con- 
straints on the supersymmetric flux. The mini- 
mum gap between the theoretical predictions and 
the data occurs at light masses and is of almost 
one order of magnitude. We have discussed by how 
much this gap changes in terms of the dark mat- 
ter distribution. This variation is relevant only for 
signals coming from the galactic center. For this 
sector, we have shown that using the r-dependent 
log-slope distribution of Ref. [l8j the gap between 
data and supersymmetric fluxes is reduced by a fac- 
tor of 3 with respect to the NEW profile. Only pro- 
files as steep as the Moore et al. one, disfavored by 
recent simulations, could exclude some light neu- 
tralino configurations. We have also shown that 
ncutralinos of masses around 30-40 GcV could ex- 
plain the EGRET excess in case of a significant 
enhancement effect as compared to the NEW dis- 
tribution. A general word of caution concerns the 
fact that the background due to conventional cos- 
mic rays production mechanisms still suffers from 
sizeable uncertainties. 

• We have shown that in the case of cosmic antipro- 
tons, no constraint on the supersymmetric param- 
eters can be derived, if one assumes a very conser- 
vative attitude in the selection of the propagation 
parameters. However, it is remarkable that indeed 
the signal at very small masses is close to the level 
of detectability. Some breakthrough in the knowl- 
edge of the astrophysical parameters could allow a 
significant exploration of small mass configurations. 
This is particularly true for neutralino masses be- 
low about 15 GeV, in view of the typical funnel 
shape displayed in the scatter plots. 

• The present measurements of up-going muons from 
the center of the Earth put some constraints on 
neutralino configurations for masses above 50 GeV. 
For lighter neutralinos, explorations by neutrino 



telescopes require a substantial increase in sensi- 
tivity with an energy threshold close to 1 GeV. In- 
vestigations of light neutralinos by up-going muons 
from the Sun are very disfavored. 

We wish to recall that, according to the measurements 
of the HEAT Collaboration the spectrum of the 

positron component of cosmic rays shows some enhance- 
ment between 7 and 20 GeV. This is only a mild effect 
which, as shown in Ref. [t^, could be explained in terms 
of conventional secondary production mechanisms. Al- 
ternatively, some authors have interpreted this effect as 
a deviation from a pure secondary flux, which could be 
generated by neutralino self-annihilation 

IzlEi- This 

hypothesis, to be a viable one, requires that the neu- 
tralino flux is enhanced by a sizable factor. Since the 
measured positrons must be created in a region around 
the Earth of a radius of a few kpc |53j . an enhancement 
would imply a significant dark matter overdensity in that 
region, with implications for the p signal. This scenario 
appears strongly model dependent, and as such not suit- 
able for setting constraints to supersymmetric parame- 
ters. 

Antideuterons in space as a signal of neutralino self- 
annihilation, which were shown in Ref. |75l| to be a very 
promising investigation means, will be considered in a 
forthcoming paper. 
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